Sub-wavelength imaging at infrared frequencies using an array of metallic nanorods 
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We demonstrate that an array of metallic nanorods enables sub- wavelength (near-field) imaging at 
infrared frequencies. Using an homogenization approach, it is theoretically proved that under certain 
conditions the incoming radiation can be transmitted by the array of nanorods over a significant 
distance with fairly low attenuation. The propagation mechanism does not involve a resonance 
of material parameters and thus the resolution is not strongly affected by material losses and has 
wide bandwidth. The sub-wavelength imaging with A/10 resolution by silver rods at 30 THz is 
demonstrated numerically using full- wave electromagnetic simulator. 

PACS numbers: 78.20.Ci, 42.30.-d, 42.70.Qs 
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I. INTRODUCTION 

Recently, the problem of sub-diffraction imaging has 
attracted a lot of interest in the scientific community. In 
2000, J. Pendry has shown that a left-handed material 
slab with negative index of refraction makes a perfect lens 
with unlimited resolution [1]. The key property that en- 
ables sub-diffraction imaging in Pendry's lens is the am- 
plification of the evanescent spatial harmonics. This phe- 
nomenon is the result of the resonant excitation of sur- 
face waves supported by the interfaces of the left-handed 
material slab. However, it soon became evident that this 
imaging mechanism is strongly sensitive to losses and dis- 
persion [2-4], besides the fact that practical implemen- 
tations of left-handed materials are still a challenge. To 
overcome these obstacles, other mechanisms to achieve 
subwavelength imaging have been studied. Namely, some 
researchers explored the fact that photonic crystals may 
enable negative refraction close to the band-gap edges 
[5, 6]. Also, recently a planar magnetoinductive lens was 
proposed for near-field imaging at the microwave range 
[7] . This approach is based on the excitation of surface 
waves at two weakly coupled resonant interfaces to am- 
plify the evanescent waves. 

A sub-diffraction near-field imaging mechanism of a 
different kind has been proposed by other researchers [8- 
14]. This technique is based neither on negative refrac- 
tion nor on the amplification of evanescent waves. The 
idea is to transport the evanescent waves through an en- 
gineered material (transmission device), rather than to 
amplify them. This canalization mechanism is possible 
if the incoming wave can couple most of its energy to 



an electromagnetic mode of the material with phase ve- 
locity nearly independent of the transverse wave vector 
(i.e. the associated wave normal surface must be nearly 
flat). In our opinion, the best solution proposed thus 
far at microwave frequencies is the one reported in [14]. 
In this configuration the engineered material is formed 
by an array of perfectly conducting wires (wire medium) 
that guide the radiation "pixel by pixel" from the in- 
put plane to the output plane. As discussed in [14, 15], 
the wire medium lens is capable of transporting subwave- 
length images without significant distortion to nearly un- 
limited distances since the influence of losses is negligibly 
small. Moreover, the resolution of the transmission de- 
vice is only limited by the spacing between the wires. The 
problem with the wire medium lens is that its design can- 
not be directly scaled to higher frequencies (infrared and 
optical domains) because metals lose their conducting 
properties, and consequently the electrodynamics of the 
equivalent medium is significantly altered. The straight- 
forward attempts to get sub- wavelength imaging at opti- 
cal frequencies using rodded medium [16] encounter prob- 
lems and limitations related to the effect of strong losses. 

Recently a non-local permittivity model was proposed 
by the authors to describe the electrodynamics of plas- 
monic rods [17]. This study enables the proper dimen- 
sioning and generalization of the transmission device pro- 
posed in [14] to higher frequencies. In [17] it was proved 
that a periodic array of rods with negative real part of 
permittivity (this is the case of all metals at infrared 
and optical frequencies) can be modeled using an effec- 
tive permittivity tensor and that spatial dispersion effects 
need to be considered. Moreover, our studies showed 
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that despite the rods not being perfectly conducting, 
the equivalent medium still supports nearly dispersion- 
less modes, and we speculated that such modes may be 
used to enable near-field subwavelength imaging at the 
infrared and optical domains. 

In this work, we will extend the theory proposed in [17] 
and demonstrate that in fact silver nanorods enable sub- 
diffraction imaging up to terahertz (THz) and infrared 
(IR) region of the electromagnetic spectrum. In recent 
years, THz and IR radiation has attracted special atten- 
tion, because it offers scientific and technological oppor- 
tunities in many fields, such as in telecommunications, 
in computer technology, in spectroscopy, and in sensing 
[18]. Near-field imaging at THz and IR frequencies may 
be very important for diverse applications, such as the 
near-field microscopy of surfaces, recording information 
with increased storage capacity, or detecting and measur- 
ing the mass and volume fractions of bio cells (e.g. cells 
of blood), many organic molecules, radioactive atoms, 
microbes, viruses and explosive matters. 

The paper is organized as follows. In section II we 
briefly review the electrodynamics of the plasmonic rod 
medium, and the possibility of canalizing waves through 
a rod medium slab. In section III we explain how to prop- 
erly design the electromagnetic crystal in order to ensure 
that the phase velocity of waves is independent of the 
transverse wave vector. Next, in section IV we introduce 
an homogenization technique that enables the accurate 
analytical characterization of the transmission properties 
of a rod medium slab. In section V it is explained how to 
tune to thickness of the transmission device in order that 
it is operated in the canalization regime. In section VI 
we present full wave simulations that illustrate the po- 
tentials of the proposed configuration, and demonstrate 
subwavelength imaging at 30THz. Finally, in section VII 
the conclusions are presented. 

In this work the fields are assumed monochromatic 
with time dependence e +Ja;t . 



II. CANALIZATION OF WAVES THROUGH 
THE ROD MEDIUM 



In order to enable the proper study and design of 
waveguiding structures based on plasmonic rods, next we 
review the homogenization theory introduced in [17]. Let 
us consider a periodic array of plasmonic rods oriented 
along the z-direction, as depicted Fig. 1 The rods are in- 
finitely long, with radius R, and the spacing between the 
rods (lattice constant) is a. The rods are embedded in a 
dielectric host with permittivity Sh- The permittivity of 
the rods is e m . 

In [17] it was proved that for long wavelengths the 
fields in the periodic structure can be characterized using 
a nonlocal permittivity dyadic: 




FIG. 1: Periodic array of long thin plasmonic rods arranged 
in a square lattice. 
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where fv = ttR 2 /a 2 is the volume fraction of the rods, 
k z is the z-component of the wave vector, (3 — uj^/EhUa 
is the wave number in the host medium, and p p is the 
plasma wave number for perfectly conducting rods given 
by [17, 19]: 



{l3 p aY 



2tt 



hi 



2ttBJ 



0.5275 



(4) 



As discussed in [17], provided the permittivity of the rods 
is such that |e m | 3> £h and the rods are thin, then e t ~ 1. 
That is the case here since the plasma frequency of metals 
is well above the far infrared band. 

It is well known [20-23], that the permittivity of met- 
als at infrared and optical frequencies can be accurately 
described by the Drude model: 



e m (w) = e 



(5) 



where u> m and T are the plasma and damping frequen- 
cies of the material, respectively. To give further insights 
into the permittivity model, let us temporarily consider 
that losses can be neglected in the frequency window of 
interest, and that u <C 0J m - In these circumstances, sub- 
stitution of (5) into (3) gives 
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e(u, k z ) = e h [e t (xx + yy) + £ zz (u>, k z )zz] 



(1) 



where we put (3 m — w my /eoHQ. By direct inspection of 
(6), it is clear that wave propagation in the rroy-plane 



3 



(i.e. k z — 0) with electric field polarized along the rods 
is only possible when (3 > /3 C ut— off 7 where: 
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(7) 



Hence, the cut-off frequency of the structure is deter- 
mined by both the plasmonic properties of the metal and 
by the plasmonic properties intrinsic to the geometry of 
the structure. 
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FIG. 2: A plane wave illuminates a slab of the rod medium 
(sideview). The thickness of the slab is L and the rods are 
normal to the interface. The structure is periodic along the x 
and y directions. 

The idea to achieve sub-diffraction imaging is to guide 
the electromagnetic fields through an array of metallic 
rods normal to the front and back interfaces, as illus- 
trated Fig. 2. The canalization phenomenon occurs 
when the phase velocity of an electromagnetic mode is 
nearly independent of the transverse wave vector. Next, 
we study in which circumstances can this happen. 

We consider that the thickness of the slab is L, and 
that for z < and z > L the medium is air. A plane 
wave with transverse wave vector k|| = (k x ,k y ,0) illumi- 
nates the slab (Fig. 2 depicts the case k y — 0, which is 
assumed throughout this section without loss of gener- 
ality). The incident magnetic field is parallel to the in- 
terface (the wave is transverse magnetic with respect to 
the ^-direction) . As is well-known, the transverse compo- 
nent of the wave vector is preserved at the boundary. The 
incident wave excites two TM-polarized electromagnetic 
modes in the rod medium with propagation constants 
k^ and k^\ respectively: 
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Note that (8) was obtained by solving the dispersion 
equation for the rod medium in the same manner as in 
[17]: it simplifies to the result (23) presented in [17] when 
t I- 

When the permittivity of the rods follows the Drudc 
model (5) and the losses are negligible, (3 C is to a first 
approximation a positive real number, (with the same 
unities as (3), that increases linearly with frequency: 
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As referred before, the canalization of waves through 
the engineered material is possible if either fc£^ or k z 2 ^ 
are nearly independent of the transverse wave vector ky 
[12]. In that case all the spatial harmonics suffer the 
same phase-shift when they propagate through the rod 
medium slab (see Fig. 2), provided the length L is also 
properly chosen [12]. In [17], two different opportunities 
to get flat isofrcquency contours were identified. It was 
shown that if either (3 C -C (3 P or (3 C » (3 P , or equivalently, 
using (10) and supposing that the permittivity of the host 
is comparable to that of free-space, if either u> <C ^ m \ffv 
or ijj » oj m \ffv, the propagation constant of one of the 
modes, given by (8), becomes nearly independent of ky. 
The two possibilities are summarized below: 

1. u -C u m VJv- In this case the nearly dispersion- 
less mode is quasi-TEM, and has characteristics 
similar to the transmission line modes of the wire 
medium [19]. The radiation is not able to effec- 
tively penetrate into the rods, and so most of the 
energy propagates in the host medium. The isofrc- 
quency contour of the pertinent mode gets flatter 
as -jj=p decreases. Note that recently [24, 25], 
it was experimentally proved that a related prop- 
agation mechanism enables waveguiding of electro- 
magnetic energy at THz frequencies. 

2. lo 3> UmVTv- In this case, the propagation mecha- 
nism is based on the surface plasmon-polariton sup- 
ported by each rod. The energy of this surface wave 
is tightly bounded to the plasmonic material, and 
consequently there is no interaction or coupling be- 
tween the rods (this is why the mode becomes dis- 
persionless). When ,? v „ increases the field con- 
finement also increases, and the isofrequency con- 
tours become more flat. 

Which of these possibilities is more interesting to achieve 
near-field sub-diffraction imaging? The answer is di- 
rectly related to the effect of losses. Indeed, in the first 
case, since the field energy is mostly concentrated in the 
host medium [25], it is expected that the rod material 
losses have a mild effect on the propagation properties. 
Quite differently, in the second case, the surface plasmon- 
polariton has a field distribution tightly bounded to the 
metallic rods [26] , and so the effect of losses may be dra- 
matic, and ultimately destroy the possibility of imaging 
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as it was demonstrated in [16]. Hence, in this work we 
will explore the first regime, for which, as referred above, 
the rods support a quasi-TEM mode. 

Firstly, we will study the influence of u m on the prop- 
erties of the modes. We consider propagation along the 
rod axes, k\\ = 0, and use the fact that et ~ 1. Then (8) 
yields, 



fc« EE fc« 



1.(2) = 1.(2) 
K, _ H z 
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Thus, the first mode has no cut-off and propagates nearly 
with the speed of light in the dielectric host, while the 
second mode can only propagate above the cut-off fre- 
quency given by (7). On the other hand, letting k\\ — > oo, 
we find that: 



^oo — K z 
1,(2) = 1.(2) 



k\ I — oo 



-JOO 



(12) 



Hence, for off-axis propagation, the attenuation of the 
second mode increases. Quite differently, the first mode 
is never cut-off and may be nearly dispersionless. Indeed, 
using (10) we can write that, 



(13) 




and so it is evident that when 

.(i) 



Pp 



is small ki 1 ^ and 



k^oo' are nearly identical, and the isofrcquency contour 
is almost flat. In this case, the mode becomes quasi- 
TEM. Notice that if the rods are perfectly conducting 
then = (3 for any k\\ and the mode is completely dis- 
persionless and exactly TEM (transmission line mode of 
the wire medium [19]). In general, to obtain good imag- 
ing at a given design frequency, the rods material shall 
be chosen such that the plasma frequency u> m is as high 
as possible. To give an idea of the possibilities, we refer 
that the plasma and damping frequencies of the follow- 
ing metals are: aluminium (Al) 3570/19.4 THz, copper 
(Cu) 1914/8.34 THz, gold (Au) 2175/6.5, and silver (Ag) 
2175/4.35 THz [20, 21, 23]. 



III. DIMENSIONING OF THE RODS 

In order to properly dimension the engineered mate- 
rial, it is also of obvious importance to assess the effect 
of the rods spacing a, and radius R, on the propagation 
characteristics. In order to study the isofrequency con- 
tours of the quasi-TEM mode, we use (4), (11) and (13) 
to obtain: 
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FIG. 3: Dispersion contours for the quasi-TEM mode at 30 
THz for silver rods with a = 215nm and various values of R. 



This result shows that a has little influence on the dis- 
persion of the quasi-TEM mode, while the radius R is 
the main factor that determines how close to unity the 
ratio k^/k^ is. Indeed, if the rods are too thin so 
that /3 m R <C 1, the isofrequency contours arc not any- 
more flat. This is illustrated in Fig. 3, where we plot 
the isofrequency contours for a medium of silver rods at 
30 THz (e Ag ~ -5143 - 746 j), for different values of R. 
The host medium is assumed to be Halcogenide glass, 
£h/£o = 2.2, and the lattice constant is a = 215nm (i.e. 
/3a = 0.2). In Fig. 3 it is seen that the contours become 
hyperbolic when the radius of the rods is very small. No- 
tice that the skin-depth in the rod material is given by 
5 rts l//3 m . Hence, the condition to get flat contours is 
equivalent to S -C R, i.e. the metal skin-depth must be 
much smaller than the radius of the wire. 

For the imaging application it is essential that not only 
the phase velocity of one of the modes is independent 
of fcii , but also that the other mode is severely attenu- 
ated so that it cannot be excited and interfere with the 

(2) 

propagating mode. Thus, it is desirable that Im{— k z '} 
is as large as possible. From (11) one can verify that 

Im{-fc^ 2) a} increases either when the lattice constant a 
decreases, or when the ratio R/a increases. As before, 
this means that very thin rods may adversely affect the 
propagation characteristics of the mode with propagation 
constant k z ' . This is illustrated in Fig. 4, where we plot 
the normalized attenuation constant as function of k\\ for 
different values of R. The rest of the parameters are as 
in the previous example. Nevertheless, the crucial factor 
that determines how severely the non-propagating mode 
is attenuated, is not the radius R but the lattice constant 
a. In fact, it can be verified that k z 2 ^ is roughly propor- 
tional to 1/a. This has a very important implication. 
Indeed, the total attenuation (in dBs) for a fixed dis- 
tance L (see Fig. 2) will be proportional to L/a. Hence, 
in general, the smaller the lattice constant is the more 
attenuated the mode is. Assuming that L is comparable 
to the wavelength in the host material, this implies that 
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FIG. 4: Normalized attenuation constant (non-propagating 
mode) at 30 THz for Silver rods with a — 215nm and different 
values of R. 



a must be chosen such that /3a -C 1. 

Let us summarize our findings. The spacing between 
the wires a must be chosen such that (3a « 1. This 
ensures that the mode in cut-off is strongly attenuated 
and does not interfere with the quasi-TEM mode. On the 
other hand, to guarantee that the phase velocity of the 
quasi-TEM mode is nearly independent of the transverse 
wave vector, we need to choose the radius of the rods, R, 
at least as large as the metal skin-depth 6. Since 2R < a, 
it is clear that the two objectives may be contradictory, 
and that a trade-off may be necessary. 

Moreover, as discussed next, the ratio R/a must also 
be kept relatively small, i.e. large values for R/a, let us 
say R/a > 0.1, cannot be admitted (also we remind the 
reader that the permittivity model (3) is only valid for 
diluted systems). The reason for this restriction is re- 
lated to the transverse impedance. Indeed, for thick rods 
the cross-sectional area occupied by the metal is large, 
and this tends to make the transverse impedance small. 
This effect is obviously undesired because it deteriorates 
the matching with the exterior, and ultimately the band- 
width of the imaging system. 

In practice, when the performance of the engineered 
material is mainly limited by the plasmonic properties of 
the rods, the best solution is to fix R slightly larger than 
the metal skin-depth, and choose a such that R/a < 
0.1. If this solution yields a relatively small value for 
(3a, then imaging with good quality can still be achieved. 
Notice that for very good conductors, the skin-depth can 
be neglected, and so it is always possible to design the 
engineered material in such a way that it has excellent 
imaging properties [14, 15]. 



IV. TRANSMISSION COEFFICIENT 

So far, the possibility of near-field imaging was dis- 
cussed mostly in a qualitative way. A more quantitative 
and detailed study requires the calculation of the trans- 



mission coefficient for the configuration shown in Fig. 2. 
As referred in [17], even though the permittivity model 
of the medium is known, the calculation of the transmis- 
sion coefficient is not trivial. The problem is that a TM- 
polarized incident wave (i.e. with magnetic field parallel 
to the interface) will excite two distinct TM-modes inside 
the rod medium slab, as discussed in section II. This phe- 
nomenon is a manifestation of strong spatial dispersion 
[27], and has an important consequence: the continu- 
ity of the tangential electric and magnetic fields at the 
interface is insufficient to calculate the scattering param- 
eters; note that for TM-incidencc, the TE-mode, with 
electric field parallel to the interface, cannot be excited 
[17]. Thus, because of the additional TM-mode in the 
artificial material slab, the system has one degree of free- 
dom. To remove this degree of freedom, an additional 
boundary condition (ABC), besides the continuity of the 
tangential electromagnetic fields, is necessary. 

In [28] an ABC was formulated and successfully vali- 
dated for the case of perfectly conducting rods. The ABC 
establishes that, 



£h,E n is continuous 



(15) 



at the interfaces, where E n is the normal component of 
the (average) electric field (i.e. E z in Fig. 2) and Eh is the 
permittivity of the host medium (which in general can 
be different in the rod medium slab and in the regions 
z < and z > L). Can such ABC be used for the case 
of non perfectly conducting rods? 

In what follows, we demonstrate that the answer is 
positive. We will show that the physical arguments used 
in [28] can be adapted to the case of plasmonic rods. Let 
us consider the structure depicted in Fig. 2. The polar- 
ization current (relative to the host medium) inside the 
plasmonic rods is given by J<j = ju; (e m — Eh) E, where 
£ m is the permittivity of the rods (for perfectly conduct- 
ing rods, the polarization current flows uniquely on the 
surface of the rods). At the interfaces the path along 
which the current flows (i.e. the rods) has an abrupt 
discontinuity. Indeed, provided |e m /eo| 3> 1 the polar- 
ization current lines are not able to penetrate into the 
regions z < and z > L, and thus the z-component 
of the polarization current must be nearly null at the 
interfaces. Moreover, as in [28], it can be proved that 
following relation holds: 



d_ . 

dz 2 



jkn ] xH av = juje h {z)E av + Ar I 3 d e^ r dxdy 



a 2 /' 

D 



( 16 ) 

where E av (z; k||) and H av (z; k||) represent the fields av- 
eraged over & z — const, cross-section of the unit cell, and 
D is the cross section of the rod. But since the tangential 
average electromagnetic fields are continuous at the in- 
terfaces and, as discussed above, the z-component of the 
polarization current vanishes at the interfaces, the above 
relation implies that the ABC (15) holds, as we wanted 
to establish (for more details see [28]). 
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Using (15), we can now easily study the transmission 
properties of the rod medium slab. Referring to Fig. 2, 
and assuming that the incident magnetic field is along 
the y-direction and has amplitude H , we can write that: 



angle of incidence is such that sin^ = fcii//3o where (3q = 
w^/£o/io is the free-space wave number. 

It can be proven that the continuity of the tangential 
electromagnetic fields and the ABC are equivalent to [28] : 



H y = Hoe 



fc<°' 



e J "-- * + pe 3k * z , z < 

~ jkxx} U2) t (2) 0<z<L 

Ter^ )(z ~ h \ z>L 

(17) 

where p is the reflection coefficient, T is the transmission 
coefficient, k^ = —j^Jk 2 — eoMo^ 2 is the longitudinal 
component of the wave vector in the free space regions, 
and A\, A2, B\, and B2 are the unknown amplitudes of 
the TM-modes in the slab. We remember that for the 
case illustrated in Fig. 2, the transverse wave vector is 
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(fez, 0,0). For a propagating incident wave, the 



[Hy] = 0, 

1 d 2 H. 



1 



dH,„ 



e t {z) dz 2 



e h (z)e t {z) dz 

+ LU 2 fJ, £h(z)Hy 



0, 

: 0. 



(18) 



In the above the rectangular brackets [...] represent the 
jump discontinuity of a given quantity at the interface 
(i.e. the quantity evaluated at the air side subtracted 
from the quantity evaluated at the rod medium side), 
and, by definition, Sh = £0 arid St — 1 at the free space 
side. Using these boundary conditions, we readily obtain 
the following linear system of equations, 
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(19) 



where wc defined ki' 1 ^ — — j*fk% 



EhUoOj 2 . In general, the solution of system (19) must be obtained numerically. 

In the particular case where the rods are embedded in air (sh = £0), and within the approximation e t = 1, the 
transmission coefficient T can be written in a compact closed analytical form (an analogous result was obtained in 
[15] for the case of perfectly conducting rods): 
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In order to validate the proposed ABC and the derived 
formulas, next we compare the results yielded by our ho- 
mogenization theory against full wave numerical simula- 
tions [29] . We consider that the rods are made of silver, 
as in section III. The frequency is 30 THz, the lattice 
constant is 215nm ((3a = 0.2), the radius is R = 21.5nm 
(R = 0.1a), and the length of the rods is L = 5.93/im 
(f3L = 1.767r). As before, the rods are embedded in a 
host medium with £h/£o — 2.2. 

In Fig. 5 and Fig. 6 we depict the amplitude and 
phase, respectively, of the transmission coefficient as a 



function of k^/(3. The star symbols represent the full 
wave simulations, whereas the dashed line represents the 
analytical results. Note that for evanescent waves the 
transmission coefficient may be larger than unity with- 
out violating the conservation of energy. As seen, the 
agreement is not very satisfactory, particularly for the 
evanescent spatial harmonics (k\\ > (3q). After extensive 
numerical simulations, we found out that apparently the 
justification is that the analytical model (3) may slightly 

underestimate the propagation constant k^P . This means 
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FIG. 5: Amplitude of the transmission coefficient as a func- 
tion of fc||//3 for silver rods at 30 THz. The "stars" represent 
the full wave result, the dashed line represents the analytical 
result, and the solid line corresponds to the analytical result 
with corrected kP . 




FIG. 6: Phase of the transmission coefficient as a function of 
kn/P for silver rods at 30 THz. The legend is as in figure 5. 



that the propagation velocity along z is slightly less than 
the velocity of light in the host medium. This second 
order effect may be caused by the effect of higher order 
multipoles. A simple way to obtain more accurate results 
is to multiply k^P by a correction factor. The correction 
factor depends, of course, on the geometry of the crystal 
and on the permittivity of the rods (for perfectly conduct- 
ing rods no correction factor is required). For the set of 
parameters of this example, we verified that the correc- 
tion factor is 1.04 (around 30 THz). The transmission 

coefficient calculated with the corrected k^ corresponds 
to the solid line in Fig. 5 and Fig. 6. As seen, the agree- 
ment becomes very good. Similar results are obtained for 
the reflection coefficient and other geometries. We stress 
that the correction factor is independent of L, and thus 
the analytical model is very convenient to make paramet- 
ric studies and understand the electromagnetic response 
of the rod medium slab. 



V. DIMENSIONING OF THE SLAB 
THICKNESS 

In this section, we will study the imaging properties of 
the system depicted in Fig. 2, and discuss how the length 
L of the rods is chosen. As is well-known, a wave at the 
front plane, z = 0, of the imaging system can be decom- 
posed into spatial harmonics. Each Fourier harmonic is 
associated with a different transverse wave vector k|| . Be- 
cause of the superimposition principle, the Fourier har- 
monics travel independently through rod medium slab. 
Thus, the transmission coefficient T can be regarded as 
a transfer function. For the ideal reconstruction of the 
image at the output plane z = L, the amplitude and 
phase of T should be ideally independent of fci | (we note 
here parenthetically that for our system this may be pos- 
sible only for the spatial harmonics with TM-z polariza- 
tion; the TE-z Fourier harmonics are not guided by the 
metallic rods). 

In order to make T independent of k\ \ , it was suggested 
in [12] to operate the transmission device under a Fabry- 
Perot resonance. More specifically, one should choose L 
in such a way that: 



k^L 



n - 



1,2, 



(21) 



Note that the integer n determines the thickness of the 
system. Provided single mode operation can be assumed 
(i.e. provided the effect of the evanescent TM-mode can 
be ignored) it is straightforward to prove that the trans- 
mission coefficient becomes T = where the "— " sign 
is chosen for n odd, and "+" sign for n even. The re- 
markable result [12] is that if k^ is independent of fcii, 
the Fabry-Perot condition is fulfilled simultaneously for 
all the spatial harmonics, and so the imaging system is 
perfect. In practice, the resolution of the system will 
be limited by the evanescent TM-mode, and also by the 
dispersion in the phase velocity of the quasi-TEM mode. 

First we analyze the effect of the evanescent mode. For 
simplicity, we will assume that the rods are embedded 
in air so that the analytical formula (20) can be used. 
Taking into account the Fabry-Perot resonance condition 
(21), and assuming that the evanescent TM-mode is very 

much attenuated across the slab, Im {-k^L} 
obtain that: 



kf ] L \ > l,wc 



f 



1 + 



_ 



(22) 



*4 



where the "— " and "+" signs are chosen as before. Since 
ki a ^ and fc^ 1 ' 2 ^ arc in general functions of k\\, the above 
expression also is, and the imaging is not perfect, even in 
the case of perfectly conducting rods where k^ is exactly 
dispersionless [15]. When the attenuation constant of 

the evanescent mode is very large, \k^\ 3> \k^\ 



,(2) | 



> | ft! 15 1, equation (22) yields T w i.e. 



and 
the 



result that we got when assuming single mode operation. 
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Let us discuss now the effect of non-flat isofrequency 
contours for the quasi-TEM mode, i.e. the consequences 

of fci 1 ^ being a function of fey. This non-ideal behavior 
means that there is no collective Fabry-Perot resonance: 
the electrical length corresponding to the resonance be- 
comes a function of ku . Thus, even if the effect of the 
evanescent TM-mode could be neglected, the imaging 
would not be perfect. Indeed, from section II we know 
that ki 1 ^ has the following bounds: 



C < 4 1] < kg 



(23) 



Hence, to recover a certain spectral component of the 
image at the output plane, L must be chosen such that it 
satisfies (21) for the pertinent kn, and consequently the 
optimum L becomes a function of h i : 



Loo < L < Lq 



(24) 



In above, Lo,oo = Rc jwr/fcoloo}- Note that L can be 

written as L = where Xh is the wavelength in the 
host material. 



2.5 



o 1.5 



I 
< 



0.5 




10 




MP 



FIG. 8: Phase of the transmission coefficient as a function of 
fe||//9 for silver rods at 30 THz. The length of the rods is L — 
qX h , where \ h = 6.74^m and q = 1.00, 0.93, 0.88, 0.87, 0.86. 
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FIG. 9: Amplitude of the reflection coefficient as a function of 
k\\/P for silver rods at 30 THz. The length of the rods is L — 
q\ h , where X h = 6.74^m and q = 1.00, 0.93, 0.88, 0.87, 0.86. 



FIG. 7: Amplitude of the transmission coefficient as a 
function of fc|i//3 for silver rods at 30 THz. The length 
of the rods is L — q\h, where Xh = 6.74pim and q — 
1.00,0.93,0.88,0.87,0.86 (each value of q is associated with 
a different curve) . 

In practice, how should one choose LI If we want to 
recuperate only the propagating waves, the choice L = Lq 
is adequate. However, in general, the condition L = Lq is 
not appropriate for the evanescent harmonics with fen > 
/?o- For these modes the Fabry-Perot condition is fulfilled 
for some length L(ku) < L . Thus, in order that near- 
held imaging is possible, we may have to make a trade- 
off. More specifically, we may have to allow some phase 
and/or amplitude distortion for the propagating waves, 
in order to improve transmission of the evanescent waves. 
Thus, the length L will be chosen in the interval defined 
by (24). 

To illustrate the discussion let us consider an example. 
Again, we consider silver rods at 30 THz, with the same 



lattice constant, host medium, and radius as in the pre- 
vious section. In Fig. 7 and Fig. 8 the amplitude and 
phase of the transmission coefficient are plotted for dif- 
ferent values of L. We used the analytical model (20) to 
plot the curves, taking also into account the correction 
factor 1.04 to better predict k^\ as discussed in the pre- 
vious section. It is seen that when the thickness of the 
rods is close to 0.93A/,, (which corresponds to the n = 2 
Fabry-Perot resonance with L = L , using the correction 
factor 1.04), the transmission coefficient is T « 1 for the 
propagating waves. However, the amplitude response of 
the evanescent waves is relatively poor. As discussed be- 
fore, the response of the evanescent spectrum can be im- 
proved by decreasing L. It is seen that when this is done 
the phase of T for the propagating waves becomes pro- 
gressively more distant from zero, even though it remains 
relatively flat. However, the amplitude response of the 
evanescent spectrum is significantly improved, and for 
L = 0.88A/J = 5.93^iin, near field imaging seems possible 
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roughly for k\\/j3 < 3.5 with relatively small phase and 
amplitude distortion. Note, that if L is greater than L 
(curve with L = l.OOA^) then the transmission charac- 
teristic has a pole very close to kn/ft f» 1. This indicates 
that for this set of parameters a surface wave is supported 
by the rod medium slab. This regime must be avoided 
since the excitation of surface waves and the selective 
amplification of very specific spatial harmonics leads to 
destruction of subwavelength imaging [15]. 

Another important aspect is to avoid interference be- 
tween the reflected wave and the source. To this end, 
we need to ensure that the reflection coefficient ampli- 
tude remains relatively small in the frequency band of 
interest. In Fig. 9 it is seen that reflection level is very 
high for L < 0.87A/j, and so the length of the rods must 
definitely be chosen larger than 0.87A^. 

To conclude this section, we study the bandwidth of 
the proposed transmission device. In Fig. 10 the trans- 
mission characteristic is depicted for different frequencies 
assuming that the length of the rods is L = 5.93/im (cor- 
responds to q = 0.88 at 30THz). It is seen that the effect 
of changing frequency for a fixed L is qualitatively similar 
to the effect of changing L for a fixed frequency (Fig. 7). 



o 



< 




KIP 



FIG. 10: Amplitude of the transmission coefficient as a func- 
tion of fc||//3 for silver rods at different frequencies. The length 
of the rods is L = 5.93/xm. 



VI. NEAR-FIELD IMAGING AT INFRARED 
FREQUENCIES 

Next, we illustrate the possibilities and potentials of 
the proposed subwavelength imaging mechanism. To this 
end, we simulated the performance of a rod medium slab 
that consists of an array of 89 x 69 silver nanorods (see 
Fig. 11) using CST Microwave Studio full-wave elec- 
tromagnetic simulator [29]. The nanorods are embed- 
ded in Halcogenide glass (e^ = 2.2e ) with total dimen- 
sions of 19.35 x 15.05 x 5.93/xm. The lattice constant is 
a = 215nm, the radius of the rods is R = 21.5nm, and the 




FIG. 11: Geometry of finite-sized slab of rod medium used in 
full- wave simulation. The critical parameters are given in nm. 



length of the rods is L = 5.93/im, consistently with the 
simulations presented in the previous sections. Namely, 
at 30THz the dispersion contours of the electromagnetic 
modes are presented in Fig. 3 and Fig. 4 (curve with 
R = 0.1a) and the transmission coefficient is shown in 
Fig. 10 for different frequencies. We stress that in all the 
simulations presented in this paper the effect of realistic 
losses in silver is taken into account. 
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FIG. 12: Distributions of \E Z \ at 29THz: a) at the source 
plane, at distance a/2 from the source and at the distance a 
from the front interface, and b) at the image plane, at distance 
a/2 from the back interface. 
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FIG. 13: Distributions of \E Z \ at 30THz: a) at the source 
plane, and b) at the image plane. 



FIG. 14: Distributions of \E X \ at 31THz: a) at the source 
plane, and b) at the image plane. 



In our numerical model the source consists of a set 
of perfectly conducting wire dipoles driven by 1A cur- 
rents in their centers and arranged in the form of let- 
ters "IR", see Fig. 11. The source is placed at the 
distance a/2 = 107. 5nm from the front interface. The 
"IR" -shaped dipole source radiates both TE-z and TM-z 
waves. As discussed before, the rod medium slab is able 
to guide only the TM-z spatial harmonics. Thus, in order 
to evaluate the waveguiding properties of the rod medium 
slab we numerically calculated the z— component of the 
electric field at the front (input) plane and at the back 
(output) plane. More specifically, the field is calculated 
at the distance a/2 from the source and back interface, 
respectively. Note, that the z-component of electrical 
field includes only the contribution of TM-z waves since 
E z = for the TE-z waves. 

In Fig. 12 \E Z \ is depicted at 29THz at the input and 
output planes, respectively. It is clear that the imag- 
ing properties are very poor, and that the letters "IR" 
are hardly perceptible in the image, even at the front 
interface. As discussed in section V, the reason for this 
phenomenon is related with the strong reflections from 
the rod medium slab due to the resonant excitation of 
surface waves. 

At higher frequencies, for example at 30THz and 
31THz, see Fig. 13 and Fig. 14, the imaging properties 
are excellent. In this frequency range the surface wave 
mode cannot be excited and so there is little interaction 



between the transmission device and the source. Due to 
this reason the "IR" shaped letters are sharply defined 
in the near-field distribution at the source plane (Fig. 
13 and Fig. 14). Also, consistently with our theoreti- 
cal model, the rod transmission device is able to canalize 
the image from the front interface to the back interface 
with subwavelength resolution. The best imaging occurs 
at 31THz. We can estimate the resolution (minimum 
distance that makes two point sources resolved) of the 
transmission device as explained next. From Fig. 10 it 
is seen that at 30THz, T > 0.7 (half-power criterion) for 
fc|| < 3.43/3, or equivalently for fc|| < 5.1(3 . Thus, the res- 
olution of the transmission device is about 5.1 times bet- 
ter than that of free-space, i.e. the resolution is around 
Ao/(2 x 5.1) = Ao/10.2. This rough criterion is consistent 
with the resolution estimated directly from the image. It 
can be checked that the thickness of the stem of letter 
"I" at the back interface is approximately 5a ss Ao/9.3. 
This implies that the resolution is around Ao/9.3. 

For frequencies larger than 31THz the imaging prop- 
erties are slightly deteriorated, as it could be anticipated 
from Fig. 10, but subwavelength resolution is still ob- 
served. This is illustrated in Fig. 15 and Fig. 16 for 
33THz and 36THz, respectively. At 36THz the image 
is not so clear as at lower frequencies, and thus we can 
treat 35 THz as the upper bound of the frequency band 
where imaging with good sub-wavelength resolution is 
observed (the lower bound is roughly at 30 THz). It is 
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thus clear that the bandwidth of the proposed waveg- 
uiding mechanism is very good: around 15%, unlike in 
other configurations proposed in the literature (e.g. [1]). 
Indeed, our conceptual approach is very different from 
other works. In our structure we are not directly explor- 
ing a resonance of the artificial material, but instead we 
are only exploiting the fact that metallic rods support a 
quasi-TEM dispersionless mode, which is by itself a very 
broadband phenomenon. 

For frequencies below 30THz the resolution and band- 
width of a properly scaled transmission device can always 
be improved [14, 15]. This is a direct consequence of sil- 
ver becoming a better electrical conductor. It is also 
worth mentioning that it is possible to get even better 
results at 30THz and higher frequencies by using alu- 
minium (Al) rods (aluminium has a larger plasma fre- 
quency than silver), or alternatively by using rods made 
of a material with an absorption peak at the frequency 
band of interest. However, it is questionable if moving to 
higher frequencies is meaningful because there are prac- 
tical limits in the realization of the nanorods. 



VII. CONCLUSION 

In this work, we assessed the possibility of near-field 
sub- wavelength imaging at infrared frequencies using ar- 
ray of thin silver rods. Using the homogenization model 



proposed in [17], we studied the dimensioning of the rods 
and of the array. It was discussed that the plasmonic 
rod medium should be operated in the regime where it 
supports a quasi-TEM mode. It was seen that subwave- 
length imaging is possible only if the isofrequency contour 
of the quasi-TEM mode is nearly flat, and, additionally, 
if the second TM-mode is strongly attenuated. In order 
to fulfill these conditions and guarantee good matching 
with free-space, it is necessary that the radius of the rods 
is larger than the skin-depth of the metal, that the lattice 
constant is much smaller than the wavelength, and that 
R/a is small. In general it may not be possible to attain 
all the conditions simultaneously, and so a compromise 
may be necessary. It was demonstrated that an addi- 
tional boundary condition is necessary to calculate the 
transmission characteristic of the rod medium slab using 
homogenization methods. It was proved that the pro- 
posed analytical formulation accurately characterizes the 
interaction of waves with the plasmonic rods, and allows 
to qualitatively understand the physics of the pertinent 
phenomena. To study the potentials and possibilities of 
the proposed subwavelength imaging mechanism, we nu- 
merically simulated the performance of an array of plas- 
monic rods excited by an "IR" shaped near field source. 
Consistently with our theoretical study, subwavelength 
imaging with resolution Ao/10 and 15% bandwidth was 
demonstrated. 
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